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Abstract. The energy finite element analysis (EFEA) is more effective for analyzing high frequency 
vibration problem compared to the traditional finite element method (FEM). When applying the 
EFEA to complicated structures, it is necessary to obtain power transfer coefficients at structural 
joints. This paper proposes the theoretical formulations of computing power transfer coefficients 
between two 3D semi-infinite solids, and demonstrates the validation of the developed theoretical 
formulations through two examples with different incident waves: dilatational wave and distortional 
wave. 


Introduction 


For high frequency vibration problem of complicated structures, the traditional finite element method 
(FEM) needs much high dense elements because the size of elements is required smaller than the 
wavelength of vibration wave. Thus, the FEM takes much long computational time for complicated 
structures. Sometimes, the FEM cannot reach accuracy results and even cannot obtain results. To 
overcome the inherent drawback of the FEM, the energy finite element analysis (EFEA) has been 
developed [1]. The primary variable--energy density in the EFEA is discontinuous at locations of 
discontinuities with a single member or at positions where different member are connected. Hence, a 
special approach based on the continuity of power flow across the joints is developed to form global 
system equation at the joints [2]. This continuity is expressed in terms of power transfer coefficients 
in the EFEA. When applying the EFEA to such structures, it is necessary to obtain power transfer 
coefficients at structural joints. 

The power transfer coefficients among plates or at the joints between plates and beam have been 
developed [3]. This paper proposes the theoretical formulations of computing power transfer 
coefficients between two 3D semi-infinite solids, and demonstrates the validation of the developed 
theoretical formulations. 


Theoretical Formulations 


For the simplicity of deriving equations, we let the x-z plane in the wave propagation plane, and y-axis 
is perpendicular to the x-z plane and outward. In the 3D solids, two types of waves: dilatational and 
distortional waves can propagate. 
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Define ® = ero as a scalar quantity that is called the dilatation, which corresponds to the 
dilatational motion of particles. Define ý =V xü as a vector quantity representing the displacement 
potential which corresponds to the rotational motion of particles. ū=wui +vj + wk is displacement 
vector. V è W = Q is satisfied for vector quantity y. 
From the partial differential equations of motion for free vibrations in a homogeneous isotropic 
elastic 3D solid, the following second-order partial differential equation (PDE) for the dilatation is 
obtained: 
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where C, = dail is wave speed of dilatational wave. 
Similarly, the following PDE for the distortional wave is obtained: 
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where C, = f is wave speed of distortional wave. 
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We suppose one dilatational plane wave (Fig. 1a) with heading to the interface @,, or one distortional 
plane wave (Fig. 1b) with heading to the interface ¢,, incidents in the upper 3D semi-infinite elastic 


isotropic solid. When the wave goes through the interface between the two 3D semi-infinite elastic 
isotropic solids, part of wave transmits in the lower 3D semi-infinite solid. Part of wave reflects back 
in the upper 3D semi-infinite solid. The transmitted and reflected waves may include both dilatational 
and distortional waves regardless the incident wave is dilatational or distortional wave. 


Fig.la A dilatational incident wave Fig.1b A distortional incident wave and 
and its transmitted and reflected waves its transmitted and reflected waves 


We assume that the scalar and vector quantities for incident, transmitted and reflected waves as 
follows: 
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(c) For the reflected wave, P, = & e „Y, = Qe 
We assume that all the components of wave numbers of waves in the x direction are the same. The 
totality of waves in the upper 3D solid is®* = $, + ®,, Y` = W, + W,.The totality of waves in 


the lower 3D solidis B = ®,,W =y,. 


There are eight unknown complex amplitudes. In order to solve the eight unknowns, we need to 
employ the six boundary conditions (displacements continuity and stresses continuity) in the interface 
(z = 0) and two distortional wave requirements [4]: Ve y* = 0,VewW = 0 „where superscript 
“+” denotes that the terms belong to the upper solid while superscript “-” means that the terms belong 
to the lower solid. 

The eight unknown complex amplitudes can be solved out by employing the eight equations: the 
relationships between strains and displacements, the relationships between stresses and strains. Once 
the complex amplitudes of waves are obtained, the displacements of particles can be calculated using 
the following equation: 


tix, y,z,t) = V®lxy,y,2z,t) + Vx Wy Yy, Z, t) (3) 
The power transfer coefficients can be calculated by the ratio of the power transmitted by the 


generated wave to the total incident wave power on the interface for a specific type wave with 
frequency @ and heading to the interface ø [3]: 
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where the indexes p , ’, i,j represent the incident wave type, the transmitted or reflected wave type, 
the medium carrying incident wave, the medium carrying transmitted wave, respectively (j=i 
represents carrying reflected wave). We use “L” to represent dilatational wave type and “S” for 
distortional wave type. 
Since the total power which is incident on the interface equals the total power which is carried away 
by the generated waves, the conservation of energy requires that >} t = 1 [3]. 


Validation of Theoretical Formulations 


The developed theory for calculation of power transfer coefficients are validated through some 
examples. 

In the examples, Poisson’s ratio isv = 0.25 . The ratio of density in the upper solid to the lower solid is 
p,/p, . The ratio of group speed of dilatational wave in the upper solid to that in the lower solid 


isC,,/C,, . Heading of the incident wave is @. 
Table 1 Power Transfer Coefficients (Dilatational Incident Wave) 


12 12 11 11 
Tir Tis TiL Tis >» T 
P/P, =1.0 Present 0.8485 | 0.0475 | 0.0454 | 0.0586 | 1.0000 
C/C =0.5 
p= 60° Ref [5] 0.8484 | 0.0475 | 0.0454 | 0.0586 | 0.9999 
L, |% 

P/P, =0.7 Present 0.9722 0 0.0278 0 1.0000 
C/C, =2.0 
p= 90° Ref [5] 0.9722 0 0.0278 0 1.0000 Lı 

L, Li 
P/P, =1.3 Present 0 0.0696 | 0.8852 | 0.0452 | 1.0000 re A 
C0 =2.0 5> 
ý =60 Ref [5] 0 0.0695 | 0.8853 | 0.0452 | 1.0000 


Table 2 Power Transfer Coefficients (Distortional Incident Wave) 


12 12 11 
Tor Tss T sr Tss b2 T 


P,/P,=1.3 | Present | 0.0951 | 0.7702 | 0.0463 | 0.0884 | 1.0000 | S,\" / 
C40, =2.0 


Ø=78.611" | Ref [5] | 0.0951 | 0.7702 | 0.0464 | 0.0884 | 1.0001 ARL, 


S; 
p,/p,=1.3 | Present 0 0.8025 0 0.1975 | 1.0000 | S4; 
C/C =2.0 ai 

i S 
ġ=90 Ref [5] 0 0.8025 0 0.1975 | 1.0000 i 
L, Si 
p,/p,=1.3 | Present 0 0.8990 | 0.0453 | 0.0557 | 1.0000 | ‘i S / 
Cafu = 2.0 s S 
9= 73.2208" | Ref[5] 0 0.8991 | 0.0452 | 0.0556 | 0.9999 


The results are listed in Table 1 for dilatational incident wave and in Table 2 for distortional incident 


12 12 12 12 ‘ 11 11 11 11 
wave. In Tables 1 and 2, 7,7, , 7,5, Ts» Ts, are transmitted power coefficients, and 7,,, Tis» Ts» Tss 


are reflected power coefficients, where the first subscript denotes incident wave type, the second 
subscript denotes the transmitted or reflected wave type. “L” represents dilatational wave type and 
“S” represents distortional wave type. 


Conclusions 


The developed theory and program can simulate wave propagation in 3D semi-infinite elastic medium 
very well. When a dilatational or distortional wave goes through the interface between two 3D 
semi-infinite mediums, both transmitted and reflected waves or part of them generate. The generated 
waves are both dilatational and distortional waves or part of them. Power carried away by the 
transmitted and reflected waves equals to the total power carried in by the incident wave as expected. 
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